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HIGHLIGHTS 


•  2D  computing  model  of  electrochemistry  coupled  with  mass  and  heat  transfer  in  DCFC. 

•  Good  agreement  between  numerical  and  experimental  results  (AAD  of  9%). 

•  DCFC  temperature  distribution  varied  by  15  °C  among  all  the  cell  sub-domains. 

•  Accurate  prediction  of  DCFC  polarization  curve  by  using  peroxide  mechanism. 

•  DCFC  output  is  sensitive  to  electrolyte  porosity  and  anode  specific  surface  area. 
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A  two-dimensional  modeling  of  a  lab-scale  planar  Direct  Carbon  Fuel  Cell  (DCFC)  of  20  mm  in  diameter  is 
developed  by  taking  into  account  of  the  electrochemical  mechanisms  and  mass  and  heat  transfer  phe¬ 
nomena  in  all  regions  of  the  cell  simultaneously.  The  electrodes  and  the  electrolyte  of  the  DCFC  are  both 
considered  as  distinct  regions  with  different  local  properties  such  as  permeability,  conductivity  and 
diffusivity.  An  improved  packed  bed  anodic  structure  with  a  finite  thickness  is  also  adopted.  General 
boundary  conditions  are  implemented  by  taking  into  consideration  the  species  concentrations  at  the 
DCFC  inlet  such  as  oxygen  concentration  which  is  a  very  important  parameter  to  determine  the  cell 
efficiency.  The  effects  of  the  main  operating  parameters  such  as  temperature,  inlet  gas  flow  velocity  and 
porosity  of  the  electrolyte  matrix  on  the  DCFC  efficiency  are  investigated.  A  sensitivity  analysis  based  on 
numerical  simulations  of  the  effects  of  cathode  kinetic  parameters  and  the  anode  specific  surface  area  is 
also  performed.  Good  agreement  is  obtained  between  numerical  results  and  experimental  data  with  an 
absolute  average  deviation  of  about  9%. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Today,  the  world  is  over-reliant  on  oil  and  natural  gas.  Solid 
fuels,  including  carbon,  provide  promising  alternative  choices 
because  of  their  wide  availability  [1],  In  addition,  developing  new 
processes  with  improved  energy  conversion  efficiency  with  the  aim 
of  reducing  C02  emissions  is  important  to  alleviate  global  warming. 
Consequently,  Direct  Carbon  Fuel  Cells  (DCFCs)  have  been  receiving 
recently  an  increasing  attention. 
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The  Direct  Carbon  Fuel  Cell  (DCFC)  differentiates  from  the 
Molten  Carbonate  Fuel  Cell  (MCFC)  and  the  Solid  Oxide  Fuel  Cell 
(SOFC)  by  using  solid  carbon  as  fuel  [2],  It  has  the  potential  to 
convert  the  chemical  energy  of  carbon  directly  into  electricity 
without  the  need  of  gasification  or  the  moving  machinery  associ¬ 
ated  with  conventional  electric  generators. 

One  of  significant  advantages  of  this  technology  is  the  high  en¬ 
ergy  density  of  solid  carbon  compared  to  that  of  alternative  fuel 
sources  such  as  hydrogen  and  methane  both  on  a  volumetric  and 
mass  density  basis.  Carbon  fuels  can  be  obtained  from  coal,  cracking 
of  hydrocarbons,  or  biochar  issued  from  biomass  pyrolysis  [3], 
DCFC  systems  show  promising  Open-Circuit  Voltages  (OCVs), 
but  the  current  densities  have  not  been  enough  to  consider  com¬ 
mercial  applications.  The  difficulty  in  carbon  oxidation  in  DCFCs  is 
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related  to  the  carbon  nature.  The  reactivity  and  fluidity  of  solid 
carbon  are  lower  than  those  of  gaseous  fuels.  Improving  the  anode 
kinetics  is  the  key  issue  to  achieve  a  high  performance  with  DCFCs 

[3] ,  Thus,  the  use  of  a  packed  bed  anode  arrangement  of  carbon 
particles  may  contribute  to  the  enhancement  of  the  anodic  kinetics 
performance  by  allowing  better  contact  between  carbon  spheres 

[4] . 

The  research  on  DCFC  technology  has  received  a  great  attention 
focusing  on  materials,  structures,  thermal  management  and  reli¬ 
ability.  In  studying  anodic  electro-kinetics,  researchers  have 
focused  on  different  aspects  of  the  anode  mechanisms  of  DCFC  and 
the  modeling  approaches  dealing  with  this  problem  are  quite 
scarce.  Few  DCFC  modeling  study  are  available  in  literature. 

Vutetakis  and  Skidmore  [5]  presented  a  complete  polarization 
model  of  a  DCFC  with  molten  carbonate  electrolyte.  The  electro¬ 
chemical  cell  employed  an  anode  of  coal  or  carbon  particles 
dispersed  in  a  molten  carbonate  electrolyte  at  500-800  °C.  Results 
are  presented  for  various  forms  and  concentrations  of  carbon 
including  open  circuit  potentials,  current  voltage  curves  and 
product  gas  evolution  rates.  The  authors  concluded  that  at  700  °C 
and  above  the  measured  potentials  were  in  agreement  with  the 
theoretical  values  of  the  C/C 02  redox  couple. 

Kornhauser  [6]  developed  a  DCFC  electrochemical  modeling 
based  on  fluidized  bed  configuration.  He  employed  molten  car¬ 
bonate,  with  carbon  fuel  particles  and  oxidizer  entrained  and 
circulated  through  the  electrodes.  One  of  the  major  achievements 
of  this  work  is  to  understand  that  the  packed  bed  geometry  is 
capable  of  producing  large  area  per  unit  volume  available  for  re¬ 
action.  Large  surface  area  when  used  for  the  chemical  reaction 
produces  high  mass  transfer  rates.  The  effects  of  temperature  and 
the  carbon  loading  on  the  performance  of  the  DCFC  are  also 
investigated. 

Liu  et  al.  [4]  presented  a  developed  ohmic  polarization  formu¬ 
lation  for  a  packed  anode  structure  of  the  DCFC.  Results  indicate 
that  anode  ohmic  polarization  is  in  the  same  magnitude  as  the 
highest  loss  affecting  the  DCFC  performance  (activation  polariza¬ 
tion)  which  is  due  to  the  low  activity  of  carbon.  The  analytical 
model  showed  also  that  approximate  total  efficiency  of  the  single 
cell  is  higher  than  55%  for  low  current  density  and  45—50%  for  high 
current  density.  Indeed,  Liu  et  al.  [7]  simulated  a  highly  efficient 
power  system  that  combines  a  methane  catalytic  decomposition 
reactor  (MCDR),  DCFC  and  SOFC.  The  methane  is  decomposed  to 
carbon  and  hydrogen.  Carbon  is  used  as  fuel  in  the  DCFC  to  generate 
power  and  pure  CO2.  Hydrogen  is  used  as  fuel  in  the  SOFC.  Elec¬ 
trochemical  modeling  of  the  DCFC  performance  was  developed  in 
order  to  predict  the  polarizations,  operating  voltage,  power  output 
and  exhaust  temperature.  This  model  presents  an  improvement  to 
the  model  described  later  [4],  Results  indicated  that  approximately 
65%  overall  cell  efficiency  was  achieved  at  a  temperature  range  of 
600—800  °C.  The  model  permitted  also  to  confirm  that  the  DCFC  is 
under  anode  activation  polarization  and  ohmic  polarization 
control. 

Li  et  al.  [8]  developed  an  anode  micro-structural  model  of  car¬ 
bon  particle  in  molten  carbonate  configuration  based  of  six  step 
mechanisms  assuming  no  rate  determining  step.  Through  this 
model,  which  accounts  of  the  irreversibility  of  the  anode  process 
and  supposing  that  the  electrochemical  oxidation  of  carbon  takes 
place  only  at  the  edges  of  the  graphene  sheets,  the  performance  of 
carbon  black  and  graphite  as  fuel  of  DCFC  was  examined  and 
showed  an  activation  loss  in  the  range  of  0.024—0.28  V  with  a 
current  density  of  10-120  mA  cm-2  at  700  to  900  °C.  The  devel¬ 
oped  analytical  model  is  based  on  a  modified  anode  reaction 
mechanism  and  a  new  method  of  rate  calculation  deduction.  This 
helps  to  clarify  and  obtain  the  kinetic  properties  of  the  DCFC  i.e. 
reaction  rate,  exchange  current  density  and  polarization  losses. 


Recently,  Cooper  et  al.  [9]  criticized  the  structure-based  kinetic 
model  of  the  anode  developed  by  Li  et  al.  [8]  and  concluded  that  it 
should  be  combined  with  porous  electrode  theory  to  study  the 
reaction  distribution  computational  models.  This  model  should  be 
also  tested  using  a  well  characterized  carbon  under  well  defined 
operating  conditions  to  ensure  the  definition  of  a  rate  equation  that 
should  be  estimated  on  the  basis  of  reasonable  determined  values 
of  the  electrochemical  factors  and  rate  constants  of  each  individual 
step  of  the  studied  mechanism. 

Chen  et  al.  [10]  proposed  a  ID  macrohomogeneous  model  of 
carbon  corrosion  in  a  DCFC  system  by  considering  4-electron  car¬ 
bon  oxidation  reaction,  2-electron  CO  oxidation  reaction  as  well  as 
the  Boudouard  reaction.  Obtained  results  show  that  the  active  zone 
is  located  in  the  interface  carbon/electrolyte  and  presents  about 
30%  of  the  total  height  of  the  carbon  particles  bed.  Simulation  re¬ 
sults  confirm  that  the  anodic  DCFC  performance  is  mainly 
controlled  by  the  electrochemical  oxidation  of  carbon. 

Apart  from  the  analytical  cited  models  on  DCFC  electrochem¬ 
istry,  a  recent  numerical  model  had  been  developed  by  Alexander 
et  al.  [11],  They  implemented  an  axisymmetric  kinetic  model  of 
carbon  bed  in  a  solid  carbon  fuel  cell  using  the  finite  element 
method  to  solve  the  coupled  kinetic,  gas  mass  conservation  and  gas 
transport  equations  inside  the  anode  bed  in  order  to  induce  the 
operational  space  of  the  fuel  cell  system.  Numerical  results  revealed 
that  the  lower  specific  surface  area  and  bulk  density  of  the  carbon 
offset  by  higher  fuel  reactivity  in  the  dry  gasification  environment 
inside  the  anode.  These  authors  concluded  also  that  the  cell  per¬ 
formance  increases  with  the  onset  of  carbon  conversion  due  to 
specific  surface  area  increase.  Then,  it  begins  to  decrease  by  around 
50%  fuel  conversion  as  the  loss  of  reactive  material  in  the  bed 
dominated  its  behavior. 

Until  now,  mathematical  modeling  of  DCFCs  is  still  at  the 
beginning  stage  and  didn’t  yet  account  for  the  coupled  electrodes 
electrochemistry  with  the  heat  and  mass  transport  mechanisms. 
Therefore,  understanding  the  interaction  between  these  complex 
phenomena  has  such  an  explorative  nature  in  this  field. 

In  the  literature,  several  concepts  of  DCFC  have  been  presented 
by  employing  conventional  electrolyte  materials  such  as  molten 
hydroxide  [1],  molten  carbonate  [12],  and  solid  oxide  [13],  The 
electrolyte  is  a  critical  material  in  the  fuel  cell.  Indeed,  the  perfor¬ 
mance  of  the  cell  is  dependent  on  the  conductivity  of  the  electro¬ 
lyte  materials.  Also,  the  electrolyte  has  to  be  dense  to  prevent  gas 
permeability  to  opposite  sides  and  has  a  large  area  to  minimize  the 
bulk  resistance.  It  should  be  an  electronic  insulator  but  a  good  ionic 
conductor.  For  state-of-the-art  DCFC  based  solid  oxide  electrolyte 
such  YSZ  needs  a  temperature  above  700  to  800  °C  to  ensure 
adequate  ionic  conductivity  in  order  to  guarantee  sufficient  power 
outputs. 

In  fact,  a  new  research  topic  has  been  devoted  to  the  develop¬ 
ment  of  new  alternative  electrolyte  materials  that  operate  at  lower 
temperatures  such  as  Samarium  Doped  Ceria  (SDC).  Ceria- 
carbonate  composite  electrolytes  have  attracted  significant  atten¬ 
tion  during  the  last  decade  for  their  application  in  Intermediate- 
Temperature  Fuel  Cells  (IT-fuel  cells). 

Composite  electrolytes  usually  show  enhanced  conductivity 
compared  to  the  conventional  solid  electrolyte  used  in  SOFCs  and 
the  molten  carbonate  electrolyte  used  in  MCFCs,  but  the  effects  are 
not  simply  additive;  rather  they  are  synergistic  which  means  that 
the  global  conductivity  is  significantly  higher  than  in  both  of  the 
constituent  phases.  The  conductivity  enhancement  is  suggested  to 
be  due  to  high  ionic  conductivity  in  the  interface  region  between 
components  [14,15],  In  addition,  composite  electrolyte  materials 
can  be  fabricated  easily  into  a  certain  shape  to  support  the  fuel  cell 
structure.  These  types  of  composite  materials  show  a  very  prom¬ 
ising  performance  when  used  as  electrolyte  in  SOFCs. 
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The  first  DCFC  apparatus  based  on  ceria-composite  electrolyte 
containing  Li/Na  carbonate  eutectic  and  doped  ceria  phase  is  tested 
byjia  et  al.  [16],  A  four-layer  pellet  cell  is  fabricated  by  a  co-pressing 
and  sintering  technique.  Commercial  activated  carbon  powder  is 
mixed  with  the  composite  electrolyte  and  is  retained  in  the  anode 
cavity  above  the  anode  current  collector.  With  a  suitable  CO2/O2 
ratio  of  the  cathode  gas,  an  operating  temperature  of  700  °C,  a 
power  density  of  98  mW  cm-2  at  a  current  density  of  200  mA  cm  2 
is  achieved. 

Elleuch  et  al.  [17]  investigated  this  type  of  composite  electrolyte 
in  a  DCFC  system  fueled  by  slurry  of  graphite  and  molten  carbonate. 
They  demonstrated  that  the  composite  electrolyte  is  promising 
compared  to  the  single  constituent’s  electrolyte. 

Although  promising  experimental  results  were  achieved  by  our 
research  group  [17,18],  design  and  operating  parameters  need  to  be 
optimized  for  maximum  performance,  so  that  a  proper  assessment 
can  be  made  for  possible  future  DCFC  technology  development.  For 
that,  a  thought  on  mathematical  modeling  and  computational 
simulations  is  recommended  in  order  to  predict  how  fuel  cell 
performance  is  affected  by  the  variation  of  operating  conditions 
and  geometrical  parameters.  Such  model  should  be  simple, 
computationally  fast  and  able  to  capture  reliably  the  DCFC  behavior. 

The  purpose  of  this  work,  and  hence  its  contribution,  is  to 
develop  a  preliminary  mathematical  model  of  a  Lab-scale  DCFC 
system  in  order  to  build  a  numerical  tool  enabling  the  numerical 
prediction  of  heat  and  mass  transport  phenomena  effect  on  the 
electrochemistry  of  the  DCFC  and  the  assessment  of  the  DCFC 
design  optimization  for  a  wide  range  of  parameters. 

This  model  takes  into  account  of  the  electrochemical,  mass  and 
heat  transfer  phenomena  in  all  regions  of  the  cell  simultaneously 
and  permits  to  describe  the  temperature  and  species  concentra¬ 
tions  profiles  all  over  the  regions  of  the  fuel  cell.  Parametric  study  is 
then  performed  to  predict  the  effect  of  key  operating  parameters 
such  as  temperature,  inlet  gas  flow  velocity,  electrolyte  matrix 
porosity  and  anodic  specific  surface  area  on  the  DCFC  performance. 
The  model  results  are  also  specifically  validated  experimentally  in 
Tianjin  Key  Laboratory  of  Applied  Catalysis  Science  and  Technology, 
China  by  using  an  actual  state-of  the-art  DCFC  system. 

2.  Experimental 

2.1.  Cell  fabrication 

Samarium  Doped  Ceria  (Ceo.sSmo.2O19)  sample  was  synthesized 
via  oxalate  co-precipitation  process.  The  SDC  precursor  is  sintered 
at  700  °C  for  2  h  to  form  the  pale  yellow  SDC  powder 
(Ceo.sSmojOl.g)  [16,17],  Then,  a  binary  carbonate  eutectic  powder, 
i.e.  Li2C03/Na2C03  in  a  mole  ratio  of  2:1,  was  prepared  by  a  process 
of  sequential  ball  milling  for  2  h,  co-sintering  at  700  °C  for  2  h  and 
grinding.  Composite  ceria-carbonate  electrolyte  material  was  ob¬ 
tained  through  mixing  the  last  two  powders  (carbonate  powder 
and  SDC  powder)  in  a  weight  ratio  of  3:7  by  ball  milling  for  2  h.  The 
composite  cathode  powder,  which  consists  of  30  wt%  the  composite 
electrolyte  and  70  wt%  the  Lithiated  Nickel  Oxide  (LixNi  i_xO)  pow¬ 
ders,  was  also  prepared  through  2  h  ball  mill  mixing  and  sintering 
at  700  °C  for  2  h.  The  cathode  has  a  porous  nature  which  facilitates 
reactant  distribution  across  the  reaction  sites  [19],  It  includes  some 
sites  rich  on  nickel.  The  main  role  is  to  accelerate  the  electro¬ 
chemical  reaction  occurrence.  More  details  the  reaction  mecha¬ 
nisms  are  available  in  Ref.  [20], 

The  two-layer  DCFC  was  fabricated  by  dry-powder-pressing 
technique  beginning  by  loading  the  composite  electrolyte  fol¬ 
lowed  by  the  composite  cathode  powder  with  the  aid  of  a  60  mesh 
sieve.  The  pressed  cell  size  was  20  mm  in  diameter  with  the  elec¬ 
trode  surface  area  of  1  cm2  and  1.4  mm  in  thickness  which  consists 


Fig.  1.  Real  design  of  the  electrolyte/cathode  pellet  used  for  the  DCFC  performance 
measurement. 


approximately  on  0.8  mm  thick  composite  electrolyte,  and  0.6  mm 
thick  composite  cathode  (Fig.  1). 

Slurry  of  anode  material  was  prepared  by  adding  molten  car¬ 
bonate  electrolyte  powders  (Li2C03/Na2C03  in  a  mole  ratio  of  2:1) 
into  graphite  carbon  particles.  Ball-milled  the  mixture  at  1000  rpm 
for  2  h,  then  the  powder-like  electrolytic  material  would  be  ready. 
The  mass  ratio  of  electrolyte  to  carbon  materials  was  chosen  to  be 
1:9  wt%.  The  carbon  particles  (graphite)  forming  the  packed  bed 
serve  both  as  fuel  and  as  electrode.  The  anode  and  cathode  current 
collectors  were  made  of  silver  ring.  The  two-layers  DCFC  was  tested 
in  a  stainless  steel  apparatus  sealed  gas  tight. 

2.2.  Experimental  setup  and  testing  procedures 

The  prepared  DCFC  pellet  shown  in  Fig.  1  was  placed  in  the 
center  of  a  tubular  electric  furnace.  A  thermocouple  (type  K)  was 
attached  to  measure  the  operating  temperature  of  the  cell.  The  cell 
was  heated  to  its  operating  temperature  at  5  °C  min-1.  Nitrogen 
was  employed  as  anode  protective  gas  and  a  mixture  of  O2  and  CO2 
was  used  as  cathode  gas  [17], 

Once  the  required  operating  temperature  was  reached,  I—V 
curves  were  collected  by  a  VersaSTAT3  Electrochemical  System 
using  VersaStudio  software  for  automatic  data  collection  based  on 
the  linear  sweep  voltammetry  method  (LSV)  at  scan  rate  of 
1  mV  s-1. 

3.  Mathematical  model  formulation 

Fig.  2  schematically  demonstrates  the  2-D  physical  domain  of 
the  DCFC  divided  into  four  sub-domains:  the  anode,  the  composite 
electrolyte,  the  composite  cathode  and  the  cathode  gas  channel. 
The  cell  size  was  1.4  mm  in  thickness  which  consists  of  a  0.8  mm 
thick  composite  electrolyte,  and  0.6  mm  thick  composite  cathode. 
The  anode  thickness  was  about  1  mm  whereas  the  cathode  gas 
channel  is  0.8  mm  thick.  The  cell  is  continuously  supplied  with 
reactants.  On  the  cathode  side,  O2  and  CO2  diffuse  from  the  cathode 
gas  channel  (bulk  phase)  to  the  cathode  reaction  site  through  the 
porous  structure.  Within  the  interface  cathode/composite  electro¬ 
lyte,  oxygen  can  reacts  alone  and  is  converted  to  oxygen  ions  (O2-) 
by  consuming  electrons  according  to  possible  Eq.  (1),  or  when 
reacting  with  carbon  dioxide,  it  is  respectively  converted  to  car¬ 
bonate  ions  (CO2-)  according  to  Eq.  (2). 
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Both  oxygen  and  carbonate  ions  could  be  transported  through 
the  dense  layer  of  ceria-carbonate  electrolyte  to  the  anodic  reaction 
site.  The  bed  is  supposed  to  be  composed  of  several  slabs  of 
spherical  carbon  particles  wetted  with  molten  salts.  Each  particle  is 
assumed  to  act  as  a  rigid  sphere  and  packed  with  a  simple  hexag¬ 
onal  pattern  [4]  (Fig.  3).  Within  the  anode  packed  bed,  carbon  re¬ 
acts  with  both  oxygen  and  carbonate  ions  currying  the  electric 
current  and  releasing  CO2  according  to  Eqs.  (3)  and  (4). 

Cathode  half  cell  reaction: 

02  +  4e“  ->  202-  (1) 

02  +  2C02  +  4e“  — >2C032"  (2) 


Anode  half  cell  reaction: 

C  +  202"  ->  C02  +  4e“ 

(3) 

C  +  2C032~  — 3C02  +  4e“ 

(4) 

Overall  cell  reaction: 


C  +  02  — ►  C02  (5) 

Without  losing  the  generic  physical  characteristics,  every  nu¬ 
merical  simulation  is  conceived  and  developed  based  on  a  set  of 
assumptions  motivated  by  a  lack  of  experimentally  evaluated 
physical  parameters.  Likewise,  in  this  model  we  assume  that  car¬ 
bonate  ions  CO2-  predominates  O2-  ions  conduction  on  the  ceria- 
carbonate  composite  electrolyte.  This  assumption  is  made  accord¬ 
ing  to  the  OCV-time  and  temperature  curve  which  indicates  the 
increase  of  the  OCV  just  after  the  melting  of  the  carbonate  phase 
[21],  So,  Eqs.  (2)  and  (4)  are  adopted  here  as  cathode  and  anode 
reaction  respectively. 


Fig.  3.  Schematic  drawing  of  the  packed  bed  anode  structure. 

As  the  cell  generates  more  current,  its  electric  potential  de¬ 
creases  due  to  the  irreversible  effects,  i.e.  the  activation  over¬ 
potentials  caused  by  the  electrochemical  reaction  energy  barriers, 
the  ohmic  overpotentials  caused  by  the  electric  resistance  of  the 
cell  components  and  the  concentration  overpotentials  caused  by 
mass  transfer  limitations. 

3.1.  Governing  equations 

There  are  different  global  methodologies  for  modeling  fuel  cells. 
The  easiest  division  to  make  is  between  macroscopic  and  micro¬ 
scopic  models.  The  microscopic  model  seeks  to  model  transport  on 
an  individual  pore  level,  whereas  the  macroscopic  ones  are  con¬ 
tinuum  and  average  over  this  level.  Although  the  microscopic 
models  may  provide  more  realistic  conditions  and  factors,  they 
require  a  lot  more  knowledge  of  the  microstructure  and  are  much 
more  expensive  in  terms  of  computation  time.  Macroscopic  models 
are  more  common  for  fuel  cells,  although  some  microscopic  details 
should  be  considered.  In  the  macroscopic  modeling  approach,  the 
exact  geometric  details  of  the  modeling  domain  are  neglected. 
Instead,  the  domain  is  treated  as  a  randomly  arranged  porous 
structure  that  can  be  described  by  a  small  number  of  variables  such 
as  the  porosity  and  surface  area  per  unit  volume.  Furthermore, 
transport  properties  within  the  domain  are  averaged  over  the 
electrode  volume.  Thus,  all  variables  are  defined  at  all  positions 
within  the  domain.  Averaging  is  performed  over  a  region  that  is 
small  compared  to  the  size  of  the  electrode  but  large  compared  to 
the  microstructure  [22], 

The  macroscopic  model  is  able  to  simulate  a  large-area  cell 
which  has  high  gas  and  carbon  utilization  and  non  uniform  tem¬ 
perature  distribution  [23],  In  fact,  this  paper  focuses  in  developing  a 
two  dimensional  macroscopic  numerical  model  based  on  the  single 
domain  approach.  For  all  regions  of  the  computational  domain 
(Fig.  2)  of  the  fuel  cell,  the  continuity,  momentum,  energy,  and 
species  transfer  equations  are  written  in  the  same  form.  However, 
the  diffusion,  convection,  and  source  terms  will  vary  depending  on 
the  local  properties  of  the  regions  (eg.  cathode  gas  channel,  com¬ 
posite  cathode,  composite  electrolyte,  and  graphitic  anode). 

The  model  is  developed  using  the  Cartesian  coordinates.  The 
pellet  thickness  is  considered  asy-direction  and  the  pellet  diameter 
as  x-direction.  The  cell  surface  area  is  equal  to  1  cm2  for  the  both 
physical  model  and  experimental  pellet. 

The  model  mathematical  formulation  assumes: 

•  Steady  state 

•  Compressible  and  laminar  flow 
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•  Isotropic  and  homogeneous  regions 

•  In  the  present  model,  the  graphitic  anode  bed  thickness  is 
assumed  to  be  constant  (1  mm).  This  assumption  is  made  on  the 
basis  of  our  experimental  observation  indicating  a  quasi¬ 
stationary  slow  drift  of  the  graphite  bed  during  the  cell  dis¬ 
charging  tests.  In  fact,  the  experimental  record  of  the  l—V  curve 
at  a  working  temperature  took  only  2  min  when  graphite  bed 
drift  velocity  is  estimated  to  about  10-5  mm  s-1  (very  slow).  The 
remaining  anode  bed  thickness  is  always  found  higher  than 
1  mm  thick. 

•  CO2  is  present  in  the  entire  computational  domain  except  in  the 
electrolyte.  In  fact,  based  on  previous  results  in  the  literature 
[4,12],  the  concentration  polarization  in  both  electrodes  caused 
by  slow  diffusion,  solution/dissolution  of  reactants  and  products 
(O2  and  CO2)  into  and  out  of  the  electrolyte,  is  neglected.  For 
that,  we  fixed  a  very  low  permeability  (10-30  m2)  for  the  com¬ 
posite  electrolyte  domain. 

•  Local  thermal  equilibrium  between  solid  structure  and  fluid 
phase  in  the  porous  medium  and  hence  one  equation  with 
effective  heat  conductivity  can  describe  the  heat  transfer  process. 

By  using  these  assumptions,  the  DCFC  mathematical  model  is 
reduced  to  the  following  transport  equations  for  the  entire  domain 
with  variable  diffusion  and  source  terms  which  will  be  described 
below  for  each  sub-domain: 

■  Conservation  of  mass 


-  =  $C02,a  +  (^02  +  f^COic 


■  Conservation  of  momentum 


djepuu)  djepuv)  3P  8  f  9(eU)\  8  /  8(eU)\  _  ep 

ax  fry  ax  ax  a*  7  ay  v  Qy  /  K 


d(epuv)  d(epw)  _  dP 

dx  +  ay  ay  ax 


■  Conservation  of  gaseous  species 


where  i  =  O2,  CO2,  JV2 


(9) 


■  Conservation  of  energy 

d(epCpTu)  djepCpTv)  8  /  8T\  3  f  dT\ 

ax  ay  i>x  v  ax  J  ay  v  ay  J  rev 


■  Conservation  of  charge 


(10) 


(11) 


3.7.1.  Cathode  gas  channel 

In  the  cathode  gas  channel,  the  continuity  Eq.  (6)  has  no  source 
terms  in  this  region  since  there  are  no  chemical  reactions  occurring. 
The  momentum  equations  (Eqs.  (7)  and  (8))  doesn’t  include  the 
darcyan  term  since  the  cathode  gas  channel  is  not  a  porous  me¬ 
dium.  The  source  terms  for  the  mass,  energy  and  potential  equa¬ 
tions  Eqs.  (9)— (11)  are  also  nil  in  this  region  of  the  computational 
domain. 

The  properties  of  the  gaseous  mixture  in  this  region  of  the 
computational  domain  are  calculated  by  taking  into  account  the 
different  species  present  locally.  The  mixture  density  ( p ),  viscosity 
(p)  and  mass  diffusion  (Di  mix)  are  respectively  given  by  Ref.  [24]: 

p  =  J2MiCi  (12) 


Epi  v— v  M 

.  >  =  E- 


(13) 


Where  />;  is  the  density,  /z,  is  the  viscosity,  M;  the  molar  mass  and  C; 
is  the  concentrations  of  the  gaseous  species  i. 


1  -y, 


Zyj/Dij 


(14) 


Where  y,-  is  the  molar  fraction  of  species  z,  given  by: 

Vi  =  q/]TCj  (15) 

The  diffusion  coefficient  in  Eq.  (14)  is  based  on  the  Stefan- 
Maxwell  equations  [24],  The  binary  diffusion  coefficient  Dy  for  two 
species  i  and  j  for  the  local  temperature  and  pressure  is  scaled  as 
given  in  Eq.  (16)  [25]: 


Dy  =  Dy(r°,p°)  (p°/p)  (r/T0)1'81  (16) 

The  specific  heat  and  thermal  conductivity  of  the  fluid  mixture 
are  given  by  Ref.  [26]: 


r  ^(QCp.O/M,- 

Cp  =  V  EQ  ’ 


i  =  02,C02,N2 


(17) 


(IS) 


Where  the  species  properties  cp 7;  and  /q  are  given  as  a  function  of 
temperature  [26], 

3.1.2.  Composite  cathode 

As  mentioned  earlier,  the  cathode  is  made  of  a  porous  composite 
material  (Lithiated  NiO-composite  electrolyte)  allowing  convective 
and  diffusive  migration  of  the  various  species.  The  continuity 
equation  Eq.  (6)  includes  mass  source  term  due  to  the  occurrence  of 
the  cathode  electrochemical  reaction. 

The  source  terms  for  oxygen  and  carbon  dioxide  species,  Eqs. 
(19)  and  (20),  account  for  the  gas  consumption  via  the  half  cathodic 
reaction. 

Oxygen  consumption: 


0o2  =  -jc/(4F) 


+  %  =  0 


(19) 
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Carbon  dioxide  consumption  in  the  cathode: 

=  -Jc/(2F)  (20) 

Chemical  reactions,  including  electrochemical  reactions,  involve 
energy  barriers  which  must  be  overcome  by  the  reacting  species. 
This  energy  barrier  is  called  the  activation  energy  and  results  in 
activation  or  charge-transfer  polarization,  which  is  due  to  the 
transfer  of  charges  between  the  electronic  and  the  ionic  conduc¬ 
tors.  Activation  polarization  and  the  cathode  local  current  density 
per  unit  volume  (jc)  is  normally  expressed  by  the  well  known 
Butler-Volmer  equation  (Eq.  (21))  as  follows  [27]: 

jc  =  clc-Joc  [exp(acr]3Ctc(neF)/RT) 

-  exp(  -  (1  -  ac)i7act,c (neF)/RT)] 


It  is  a  function  of  the  mixture  diffusion  of  the  gaseous  species  in 
all  regions  of  the  computational  domain  which  is  calculated  by 
taking  into  account  the  different  species  present  locally  Eq.  (25). 

D,.mix  =  (i  -y,)/J2yj/Dv  (25) 

i*i 

The  energy  equation  is  expressed  by  Eq.  (10)  and  contains 
source  terms  for  reversible  heat  Eq.  (26),  heat  produced  via  acti¬ 
vation  losses  Eq.  (27),  and  ohmic  heat  resistance  generated  due  to 
the  irreversible  resistance  (material  resistance)  to  current  flow  Eq. 

(28). 

<W  =  Jc(TAScath0de/4F)  (26) 

^act,c  =  jc'1?act,c  (27) 


ac  is  the  transfer  coefficient  considered  to  be  the  fraction  of  the 
change  in  polarization  that  leads  to  a  change  in  the  cathodic  reac¬ 
tion  rate  constant  [27-30],  Meanwhile,  the  used  cathode  charge 
transfer  coefficient  (ac)  values  in  the  present  model  are  between 
0.75  and  0.25  dependently  on  the  applied  potential  (Vceii).  The 
variation  of  the  cathodic  charge  transfer  coefficient  is  proposed  due 
to  the  possible  non  uniform  current  or  potential  distributions  in 
porous  composite  cathode  structure  which  may  be  due  to  the 
presence  of  semiconducting  layers  at  the  cathode/electrolyte 
boundary.  The  heterogeneity  of  the  composite  cathode  is  seen 
clearly  in  the  SEM  micrographs  [17]  reflecting  a  random  mixture  of 
pore  diameters  and  lengths  resulting  in  a  difference  in  the  con¬ 
ductivity  from  one  region  to  the  other.  The  value  of  ac  less  than  0.5 
are  seen  also  for  the  oxygen  reduction  reaction  through  a  com¬ 
posite  cathode  (LSM— YSZ)  in  SOFC  systems  at  high  overpotentials 
[28-30], 

The  cathode  active  surface  area  ( a_c )  is  also  known  as  the  spe¬ 
cific  reactive  surface  area  (m2  m  3).  j0c  is  the  exchange  current 
density  of  the  cathode.  It  is  detailed  in  Eq.  (22)  based  on  the 
peroxide  mechanisms  [4], 

Joe  =  ig.cpco2xPo2  (22) 

The  momentum  Eqs.  (7)  and  (8)  are  similar  to  those  developed 
in  the  cathode  gas  channel,  but  contain  an  additional  source  term 
which  is  based  on  the  Darcy  law  traducing  the  pressure  drop  caused 
by  frictional  drag  in  the  porous  cathode.  This  term  is  function  of  the 
medium  permeability  (Kcathode).  The  permeability  given  by  Eq.  (23) 
is  defined  as  a  function  of  the  porosity  (eCathode)  and  the  pore 
diameter  (Dp)  of  the  porous  cathode  medium  [31], 

^cathode  =  (jJpecathode)  j ^150(1  —  t'cathode)  )  (23) 


The  pore  diameter  is  normally  between  8  pm  and  12  pm  for 
MCFC  electrodes.  Thus,  the  permeability  is  estimated  to  be  around 
1.9  10-12  m2  [32], 

The  species  conservation  Eq.  (9)  and  the  energy  equation  Eq. 
(10)  keep  the  same  form  as  those  given  in  the  cathode  gas  channel 
compartment  but  with  up-scaled  variables  (velocities  and  con¬ 
centrations)  and  corrected  fluid  property  (diffusion  coefficient)  to 
account  for  the  porous  domain. 

Where  Djmix.eff  is  the  effective  mixture  diffusivity  coefficient 
taking  account  of  the  solid  phase  fluid  presence  and  given  by  Eq. 
(24). 


mi 


(24) 


^ohm,c  =  JcMat  (28) 

For  the  reversible  heat  source  term,  the  entropy  change  related 
to  the  half  cathodic  reaction  is  taken  from  the  MCFC  literature  due 
to  the  same  reaction  occurring  in  the  cathode  electrode  (Eq.  (2)) 
[33], 

The  effective  thermal  conductivity  of  the  porous  cathode  side  is 
given  by  Eq.  (29). 

^eff, cathode  **  ecathode^f  +  (J  —  ecathode)^s  (29) 

Where  7f  is  the  thermal  conductivity  of  the  fluid  mixture  given  by 
Eq.  (18)  and  Xs  is  the  solid  thermal  conductivity  of  the  composite 
cathode  material. 

The  potential  conservation  equation  Eq.  (11)  include  a  source 
term  given  by  Eq.  (30)  representing  the  consumption  of  ions  via  the 
electrochemical  reactions  in  the  cathode  side.  It  depends  on  the 
current  density  which  is  defined  using  the  Butler-Volmer  relation, 
Eq.  (21). 

%  =  2F0Ol  =  jc/2  (30) 

The  mass  and  heat  transfer  balances  depends  on  the  cathodic 
reaction  rate  which  is  based  on  the  value  of  the  local  electrical 
current  density.  It  is  evaluated  by  the  microscopic  approach  based 
on  the  peroxide  mechanism  in  order  to  describe  the  electrochem¬ 
istry  of  the  composite  cathode  [4,23,32], 

3.1.3.  Composite  electrolyte 

The  electrolyte  is  a  mixture  of  molten  carbonate  with  SDC 
powder  (30/70  wt%  =  Li2Na2C03/SDC)  sandwiched  between  the 
two  fuel  cell  electrodes.  The  SDC  powder  is  usually  used  as  a  porous 
solid  support  matrix  to  impregnate  the  melting  carbonates.  It  is 
capable  of  transferring  ions  from  the  cathode  to  the  anode  and  it 
blocks  the  gaseous  species  transfer  such  as  02,  N2,  and  C02. 
Therefore,  the  permeability  for  the  Darcy  term  is  set  to  be  a  very 
small  number  ( K  —  10-30  m2)  in  order  to  inhibit  the  migration  of  the 
gas  species  due  to  pressure  drop  between  the  cathode  and  the 
anode  of  the  DCFC.  Since  the  approach  chosen  here  is  the  single 
domain,  the  momentum  equations  are  solved  in  the  electrolyte  and 
they  result  into  negligible  velocity  across  it. 

Eqs.  (10)  and  (11 )  describe  respectively  the  energy  conservation 
and  the  ionic  transport  in  the  electrolyte.  The  ionic  source  term  in 
the  electrolyte  is  nil  since  no  chemical  reaction  is  occurring. 

The  energy  equation  in  the  electrolyte  is  similar  to  the  one 
developed  in  the  cathode  side,  but  account  for  the  local  thermal 
conductivity  change,  the  electrolyte  porosity  and  the  ohmic  heat 
dissipation  due  to  the  ionic  resistance.  The  ohmic  source  term  is 
given  by  equation  Eq.  (31). 
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1000/T  (K-1) 


Fig.  4.  Conductivity  of  the  mixture  of  67:33  mol%  Li2C03-Na2C03/SDC  in  3/7  weight 


‘/’ohm.elec  —  JP/ffelec  C31) 

The  average  direct  carbon  fuel  cell  current  density  is  calculated 
using  the  local  ionic  current  through  the  electrolyte  and  is  calcu¬ 
lated  using  the  ohm’s  law  at  the  interface  cathode/electrolyte: 


The  conductivity  of  the  composite  electrolyte  is  deduced  from 
the  EIS  data  and  written  as  follows  [34]: 


°eiec  =  A/Texp(-Ea/RT)  (33) 

Temperature  dependence  of  the  conductivity  for  the  30/70  wt 
%  =  Li2Na2CC>3/SDC  powders  is  plotted  in  an  Arrhenius  relation  as 
shown  in  Fig.  4. 

The  measured  conductivity  of  the  composite  electrolyte  is  be¬ 
tween  10-3  and  10-1  S  cm  1  in  the  temperature  range  of  350— 
700  °C  much  higher  than  that  of  pure  SDC.  The  activation  energy 
(£a)  of  the  composite  electrolyte  is  derived  by  plotting  the  log  (oT) 
as  function  of  (1000/T)  using  Arrhenius  expression.  The  obtained 
activation  energy  is  dependent  of  temperature  range  and  is  low  in 
both  ranges  supporting  the  high  ionic  conduction  and  trans¬ 
portation  of  the  used  composite  electrolyte  (Table  1 ). 

3.1.4.  Anode 

The  anode  consists  of  carbon  particles  (graphite)  that  serves 
both  as  fuel  and  as  electrode.  Some  carbon  particles  that  are  in 
contact  with  the  electrolyte  forms  a  paste  of  carbon  fuel  granules  in 
a  molten  salt.  In  this  region  of  the  cell,  the  carbon  electro-oxidation 
reaction  takes  place  generating  heat. 

The  governing  equations  are  similar  to  those  developed  in  the 
cathode  side.  The  continuity  source  term  Eq.  (34)  is  taken  into  ac¬ 
count  due  to  the  production  of  CO2  in  the  anode  compartment. 


Table  1 

Measured  conductivity  of  30/70  wt%  =  (67:33  moU  Li2C03-Na2C03/SDC). 
Temperature  region  Activation  energy  (£a)  Pre-exponential  factor  (A) 
(°C)  (I<J  mol-1)  (eV)  (S  cm"1  K) 

350-500  12.49  0.13  9.79E3 

500-700  38.6  0.4  1.1 2E7 


<2>rev,a  =  Ja(3TASanode/4F)  (36) 

<Z>act,a  =  ja-1?act,a  (37) 

The  reversible  heat  source  term  in  the  reaction  sites  of  the 
electrolytic  graphite  bed  is  dependent  on  the  entropy  change 
associated  with  the  DCFC  anodic  reaction.  This  term  cannot  be 
directly  calculated  by  thermodynamic  correlations,  because  they 
involve  the  entropy  associated  to  the  ionic  and  electronic  phase 
[35],  Certain  authors  tried  to  determine  one  of  the  entropy  changes 
experimentally  and  to  calculate  the  second  as  the  difference  be¬ 
tween  the  entropy  associated  to  the  overall  reaction  (Eq.  (5))  and 
the  known  entropy  (in  this  case,  the  cathodic  reaction  entropy 
AScathode).  For  the  DCFC  based  molten  carbonate  electrolyte,  the 
same  cathodic  reaction  is  employed.  So,  a  literature  research  leads 
to  the  value  of  AScathode  —  -216.2  J  mol-1  K-1  [33],  The  overall 
reaction  entropy  change  of  the  DCFC  system  is  about 
ASoveraii  =  1.6  J  mol-1  K-1  at  ambient  temperature.  In  fact,  after 
taking  into  account  of  the  stoichiometric  coefficient  we  can  notice 
that  the  entropy  change  of  the  anodic  reaction  within  the  DCFC  is 
around  ASanode  —  217.323  J  mol-1  K_1. 

The  dissipation  of  heat  in  the  DCFC  is  governed  not  only  by 
electrochemical  losses  but  also  by  reversible  losses.  Most  of  the 
reversible  heat  is  dissipated  in  the  cathode  side.  On  the  other  hand, 
anode  reaction  has  a  positive  entropy  change  leading  to  the  con¬ 
version  of  heat  into  power. 

The  potential  conservation  equation  Eq.  (11)  include  a  source 
term  given  by  Eq.  (38)  representing  the  production  of  ions  via  the 
electrochemical  reactions  in  the  packed  bed  anode. 

%  =  (2F/3)<Z>co2,a  (38) 

It  is  dependent  on  current  density  which  is  defined  using  the 
Butler— Volmer  relation,  Eq.  (39). 

ja  =  a-Q-ioa  [exp(aaf)act:a(neF)/FT) 

-  exp(  -  (1  -  a3)r,3CtA(neF)/RT)]  ' ' j 

The  activation  loss  is  dependent  on  charge  transfer  coefficients 
(aa)  and  on  the  anode  active  surface  area  (a_a).  The  used  anodic 
charge  transfer  coefficient  (aa)  was  0.75. 

The  anode  exchange  current  density  (/oa)  is  related  to  the  charge 
transfer  resistance  Rct  (Eq.  (40))  and  is  given  by  Eq.  (41)  [4], 

Ret  =  30.79  +  1 .48  exp(T/28.58)  (40) 


<Z>co2,a  =  3ja/(4f)  (34) 

The  energy  equation  depends  on  the  ohmic  resistance  Eq.  (35), 
reversible  heat  Eq.  (36)  and  heat  generated  from  the  activation  loss 
Eq.  (37). 

<2>ohm,a  =  Ja/ °c  (35) 


joa  =  RT/4FRct  (41) 

One  of  the  objectives  of  the  present  work,  unlike  previous 
studies,  is  to  provide  all  required  model  input  parameters  to  help 
scientists  to  have  a  better  understanding  of  the  model  formula¬ 
tions.  Hence,  Tables  2  and  3  exhibit  all  physical  and  thermal,  kinetic 
properties  of  the  materials  involved  in  the  mathematical 
formulation. 
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Table  2 

Operating  conditions  [17], 


Description 


Anode  inlet  pressure  (atm)  1 

Cathode  inlet  pressure  (atm)  1 

Fuel  cell  temperature  (K)  973 

Length  of  the  DCFC  (m)  IE-2 

Width  of  the  DCFC(m)  IE-2 

Carbon  dioxide  inlet  flux  (ml  min-1)  120 

Oxygen  inlet  flux  (ml  min-1)  60 

Cathode  inlet  flow  velocity  (m  s-1)  0.31 

Cathode  gas  channel  thickness  (m)  IE-3 

Cathode  thickness  (m)  0.6E-3 

Electrolyte  thickness  (m)  0.8E-3 

Anode  thickness  (m)  IE-3 


3.1.5.  Boundary  conditions 

The  applied  boundary  conditions  are  as  follows: 

x-direction 

•  For  the  cathode  inlet  (x  =  0, 0  <  y  <  ygc),  a  parabolic  profile  for  u 
with  average  cathode  gas  velocity  uc  is  specified.  The  gaseous 
species  flow  the  cathode  gas  channel  only  following  the  x-di- 
rection.  In  the  inlet  boundary  of  the  cathode  gas  channel  where 
only  O2,  N2  and  CO2  are  flowing,  there  is  no  ion  flow.  The  tem¬ 
perature  of  the  gaseous  species  in  the  inlet  of  the  cathode  is 
imposed  as  Tc.  These  conditions  are  described  by  the  following 
set  of  equations: 

v  =  0 

Qx)2,a  =  0;  Qd2,c  =  Cl!  Co2  =  C2,  Cn2  =  C3 

dcp/dX  =  0  1  ’ 

T  =  Tc 


For  the  cathode  side,  electrolyte  and  anode  side  inlet  boundary 
conditions  ie.  at  {x  —  0,  ygc  <  y  <  Y),  we  suppose  an  insulated  and 
impermeable  boundary.  Also  no  ion  flow  is  present  in  this  inlet 
boundary: 


u  =*  y  c=  0 

dT/dx  =  0,  dCi/dx  =  0  and  3 cp/dx  =  0 


(43) 


Table  3 

Parameters  used  in  the  simulations. 


Symbol 

Parameter 

Value 

References 

Panode 

Anode  inlet  pressure  (atm) 

1 

[4,18] 

P cathode 

Cathode  inlet  pressure  (atm) 

1 

[4,18] 

Tc 

Fuel  cell  inlet  temperature  (I<) 

973 

[4,18] 

Vth 

Open  circuit  voltage  (V) 

1.02 

[44] 

Y 

Length  of  the  DCFC  (m) 

IE-2 

[4] 

Z 

Width  of  the  DCFC  (m) 

IE-2 

[4] 

Pc  02,c 

C02  partial  pressure  in  the  cathode 
inlet  (atm) 

0.333 

[17] 

po2 

Oxygen  partial  pressure  in  the  cathode 
inlet  (atm) 

0.666 

[17] 

PC02  ,a 

Produced  anodic  C02  partial 

1  atm 

pressure  (atm) 

Anode  charge  transfer  coefficient 

0.75 

Cathode  charge  transfer  coefficient 

0.25-0.75 

W-* *o 

Thermal  conductivity  of  the 
lithiated  NiO 
cathode  (W  m-1  K-1) 

0.9 

[32] 

As,SDC 

Thermal  conductivity  of  the  SDC 
(W  m-1  K-1) 

1.2 

[45] 

o  Porosity  of  the  lithiated  NiO  cathode 

0.62 

[32] 

esoc 

Porosity  of  the  SDC/carbonate 
composite  electrolyte 

0.3524 

[46] 

Dp 

Pore  diameter  of  the  lithiated 

NiO  cathode  (pm) 

8-12 

[32] 

Xrathode 

Permeability  of  the  cathode  (m2) 

1.6E-12 

[32] 

Specific  area  of  the  cathode  (m-1) 

3E5 

[32] 

£c 

Concentration-independent  exchange 

500 

[4] 

current  density  (Am-2) 

Power  of  C02  partial  pressure 

-1.25 

[4] 

Power  of  O2  partial  pressure 

0.375 

[4] 

Conductivity  of  the  lithiated  NiO 
cathode  (S  m-1) 

3800 

[47] 

W„lyte 

Thermal  conductivity  of  the  electrolyte 
(W  m-1  K-1) 

2 

[32] 

w 

Thermal  conductivity  of  the  anode 
(W  m-1  K-1) 

8-12 

Porosity  of  the  packed  bed  anode 

0.396 

HI 

Specific  area  of  the  anode  (m-1) 

6E3 

[10] 

<Tc 

Conductivity  of  carbon  particles  (S  m-1) 

1.6E5 

[4] 

flow.  For  the  use  of  the  Butler-Volmer  model,  the  electronic 
potential  in  the  cathode  is  supposed  nil: 

U  ==■  v  =  0 

3T/0y  =  0,  0Q/ay  =  0  and  dcp/dy  =  0  (46) 

q>ec  =  0 


For  the  exhaust  cathode  gas  channel  boundary  conditions  at 
(x  =  X,  0  <  y  <  ygc),  there  are  a  fully  developed  flow  with  no  ion 


From  the  upper  boundary  conditions  i.e.  at  the  anode,  we  sup¬ 
pose  that  the  CO2  is  realized  following  a  fully  developed  flow. 
For  the  use  of  the  Butler-Volmer  model,  the  electronic  potential 
in  the  anode  is  supposed  as  follow: 


du/dx  =  0  and  dv/dx  =  0 

dT/dx  =  0,  dQ/dx  =  0  and  dcp/dx  =  0 


(44) 


For  the  exhaust  boundary  conditions  of  the  anode  side,  elec¬ 
trolyte  and  cathode  meaning  at  (x  =  X,  ygc  <  y  <  Y),  we  suppose 
the  same  boundary  conditions  for  the  inlet: 


dT/dx  =  0,  dQ/dx  =  0  and  dcp/dx  =  0  1 

y-direction 

•  For  the  lower  boundary  conditions  i.e.  at  the  cathode  gas 
channel,  the  following  boundary  conditions  are  applied  showing 
the  use  of  an  impermeable  and  insulated  boundary  with  no  ion 


3u/3y  =  0  and  dv/dy  =  0 

dT/dy  =  0  and  dcp/dy  =  0  .  . 

Cco2,u  =  C4,  CC02,c  =  0,  C02  =  0,  CN2  =  0  l4/j 

</>e,a  =  '/th  -  '■''cell 
3.1.6.  Numerical  approach 

The  set  of  differential  equations  governing  the  DCFC  model 
along  with  the  boundary  conditions  given  above  are  solved  by  us¬ 
ing  the  Finite  Volume  Method  based  on  the  SIMPLE  algorithm  [36]. 
The  equations  are  solved  using  a  non  uniform  staggered  grid.  The 
number  of  grid  points  used  in  the  x-direction  and  y-direction  are  22 
by  45  respectively.  The  mesh  is  considered  to  be  fine  enough  to 
capture  all  the  transport,  chemical  and  electrochemical  processes 
occurring  in  the  cell  especially  in  the  complex  cathode  Three- 
Phase-Boundaries  (TPB)  and  in  the  anode  bed.  Cell  component 
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dimensions  can  be  found  in  Table  2.  The  FORTRAN  code  in  the 
present  model  uses  the  power  law  scheme  based  on  the  Peclet 
number  for  calculating  the  convective  and  the  diffusive  flux  in  each 
differential  equation.  In  principle  the  model  is  developed  for  3D 
problems.  However,  since  only  one  mesh  element  is  generated  in 
the  z-direction,  it  becomes  2D  in  this  study. 

There  are  generally  two  ways  of  solving  for  the  electrochemical 
cell  efficiency:  either  the  operating  current  density  is  given  and 
different  potential  losses  are  calculated  or  the  cell  potential  is  set 
and  the  current  density  delivered  by  the  fuel  cell  is  then  calculated 
using  Eq.  (32).  The  latter  is  called  voltage  to  current  approach  (VTC) 
which  has  been  adopted  in  this  study  [32], 

4.  Results  and  discussions 

4.1.  Model  validation 

The  model  results  are  validated  with  the  experimental  results 
carried  out  for  the  design  and  operating  conditions  given  in  Table  2 
and  the  simulation  parameters  given  in  Table  3. 

The  experimental  measurements  of  the  DCFC  performance  were 
performed  at  700  °C  using  a  cathode  gas  mixture  of  120  ml  min  1 
CO2  flux  and  60  ml  min-1  O2  flux.  The  N2  flux  is  used  as  gas  carrier 
in  the  present  model.  The  same  electrode  material,  thickness  are 
employed  in  the  modeling.  In  addition,  the  adequate  composite 
electrolyte  conductivity  was  estimated  using  Electrochemical 
Impedance  Spectroscopy  (EIS)  method  due  to  the  large  variations 
observed  for  several  Ceria-carbonate  composite  electrolyte 
samples. 

As  shown  in  Fig.  5,  the  modeling  results  are  in  good  agreement 
with  the  experimental  measurements  over  the  entire  current 
density  range  with  an  absolute  average  deviation  (AAD)  of  about 
9%.  The  consistency  of  these  two  curves  by  itself  is  a  model- 
validation  measure  but  it  also  predicts  a  modeling  of  the  electro¬ 
chemical,  mass  and  heat  transport  on  all  the  DCFC  sub-domains.  For 
the  high  current  density  range  the  present  model  gives  accurate 
prediction  of  the  real  cell  behavior  (AAD  within  3%)  and  at  low 
current  density  range  the  AAD  is  within  15%.  These  errors  may  be 
due  to  the  fact  that  the  activation  polarization  is  modeled  using  an 
exponential  law  (Butler— Volmer  equation)  whereas  during  the 
experiments  a  linear  behavior  at  high  voltage  and  low  current 
densities  is  depicted.  It  seems  like  the  linear  ohmic  losses  is 
majoring  the  DCFC  polarization  curve  over  the  activation  and 
concentration  losses.  A  detailed  polarization  analysis  will  be 
developed  in  the  next  section. 

For  the  following  discussions,  the  parameters  given  in  Tables  1 
and  2  are  used  as  the  default  values  and  a  parametric  study  will 
be  performed. 

4.2.  Simulation  results 

Fig.  6  shows  the  different  overpotential  losses  (77)  affecting  the 
performance  of  the  fuel  cell  at  700  °C  as  function  of  the  current 
density.  The  polarization  curve  of  the  DCFC  seen  in  Fig.  5  seems  to 
have  a  linear  behavior.  This  is  suggested  to  be  related  to  the 
domination  of  the  ohmic  polarization  [17],  Conversely,  opposite 
behavior  is  observed  in  the  simulated  overpotentials  shown  in  the 
same  figure.  At  highest  current  densities  (i.e.,  higher  Vth  -  \Zceii ),  the 
ohmic  overpotential  become  constant  and  the  increase  in  losses  is 
solely  due  to  the  anodic  and  cathodic  activation  overpotentials 
which  are  directly  related  to  the  electrochemical  reaction  kinetics. 
At  high  cell  voltage  (Vceii),  the  dominant  overpotential  is  the  anodic 
activation  polarization.  As  the  cell  voltage  drops,  the  current  pro¬ 
duced  by  the  DCFC  increases,  the  cathodic  activation  polarization 
become  more  important.  The  cathode  activation  loss  exhibits  a 


linear  function  of  cell  voltage  since  the  current  density  is  a  hyper¬ 
bolic  function  of  the  cell  voltage  as  indicated  by  Fig.  5.  The  com¬ 
bined  effect  of  both  electrode  activation  loss  results  to  the  linear 
polarization  curve. 

As  it  is  shown  the  anode  and  cathode  concentration  over¬ 
potential  doesn’t  affect  the  cell  voltage  curve  at  high  cell  voltage 
drop  due  to  the  fact  that  the  anode  is  fueled  by  solid  carbon  which 
is  always  present  during  the  cell  test  period.  The  issue  was  found  to 
be  in  the  porous  composite  cathode  where  the  availability  of  gas 
mixture  is  always  crucial  which  should  neglect  the  drop  off  in  the 
amount  of  current  generation  by  the  mass  transfer  limitation 
process. 

The  temperature  distribution  along  the  DCFC  shown  in  Fig.  7  is 
simulated  for  VCeii  =  0.4  V  corresponding  to  the  voltage  giving  the 
highest  power  density  in  this  case.  On  the  light  of  Fig.  7,  it  is  evident 
that  the  temperature  gradient  in  the  x-direction  is  more  significant 
compared  to  that  in  the  y-direction  due  to  the  small  dimensions  in 
the  y-direction  compared  to  the  x-direction.  It  can  be  seen  that  the 
heat  of  the  gas  is  transferred  along  the  x-direction  via  conduction 
and  along  the  y-direction  mainly  by  convection  of  gas  flow. 

In  addition  to  the  electrolytic  ohmic  loss,  the  major  energy  loss 
in  the  DCFC  is  due  to  electrode  surfaces.  There  is  a  rapid  reduction 
observed  which  has  been  attributed  to  the  activation  over¬ 
potentials  of  both  electrode  domains.  Note  that  most  of  the  heat 
dissipated  in  this  case  is  due  to  the  activation.  Since  activation 
magnitude  is  similar  in  the  anode  to  the  cathode  (0.24  V)  at  this 
cell  voltage,  the  heat  is  generated  similarly  within  the  anode  and 
cathode  regions. 

The  temperature  in  the  cathode  region  is  slightly  lower  than  the 
anode  due  to  the  reversible  heat  which  is  converted  into  electrical 
power  in  the  anode  side  by  assuming  a  positive  entropy  change  of 
this  electrochemical  reaction  site  (ASanode  >  0). 

In  general  the  temperature  variation  (AT)  from  the  inlet  973 1<  to 
the  outlet  987.6  K  seems  to  be  around  15  °C  which  could  be 
explained  by  the  near  zero  entropy  change  of  the  single  DCFC 
system. 

Fig.  8  shows  the  carbon  dioxide  (CO2)  profile  yield  in  the  anode 
sub-domain  of  the  DCFC.  The  obtained  results  show  that  the  high 
carbon  dioxide  concentration  value  is  present  in  the  interface 
anode/electrolyte  where  the  carbon  anodic  oxidation  to  carbon 
dioxide  occurs.  The  concentration  of  CO2  diffusing  within  the 
porous  anode  decreases  slightly  through  the  anode  packed  bed 
thickness.  It  is  obvious  that  the  upper  boundary  condition  for 
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Vth-Vcell(V) 

Fig.  6.  DCFC  potential  losses  at  700  °C. 


carbon  dioxide  concentration  is  set  at  atmospheric  pressure  and 
anodic  fuel  cell  temperature  to  13  mol  m-3. 

Fig.  9  shows  the  oxygen  concentration  (O2)  profile  in  the  cath¬ 
ode  sub-domain.  The  oxygen  concentration  starts  from  the  inlet 
value  and  decreases  through  the  gas  channel  and  the  porous 
composite  cathode  media  as  it  is  consumed.  The  effect  of  the 
consumption  of  oxygen  can  be  seen  in  the  three  phase  boundary 
(TPB)  of  the  cathode  (the  higher  concentration  gradient).  This  could 
be  explained  by  the  local  cathode  electrochemical  reaction  occur¬ 
rence.  It  is  also  notable  that  the  oxygen  concentration  gradient  is 
significant  in  both  directions  due  to  the  O2  2D-diffusion  and  the 
continuous  consumption  in  the  reaction  zone. 

4.3.  Parametric  study 

In  this  section,  we  propose  to  investigate  the  effect  of  some  key 
parameters  (temperature,  anode  specific  surface  area,  cathode  ki¬ 
netics,  cathode  gas  composition  and  composite  electrolyte 
porosity)  on  the  DCFC  performance. 

The  validated  model  was  used  to  perform  in  depth  analyses  to 
gain  understanding  of  the  activation,  concentration  and  ohmic 
overpotentials  of  a  DCFC  under  different  operating  conditions. 

Fig.  10  shows  the  effect  of  temperature  on  DCFC  performance. 
Typically,  the  operating  temperature  of  the  DCFC  is  above  the 
melting  temperature  of  the  carbonate  in  the  range  of  600-700  °C. 
The  modeling  results  predict  well  the  experimental  data  depen- 
dently  on  the  operating  temperature.  As  the  operating  temperature 
increases,  an  improvement  in  the  fuel  cell  performance  is  observed 
due  to  the  fact  that  physico-electrochemical  properties  such  as 
activation  energies,  binary  species  diffusivities  and  ionic  conduc¬ 
tivities  are  enhanced  with  the  increase  of  temperature.  This  effect  is 
clearly  illustrated  in  Fig.  10  by  the  slope  of  the  polarization  curve 
that  decreases  from  364  pQ  m-2  at  873  I<  to  269  pQ  m-2  at  973  K 
more  at  higher  temperatures. 

The  activation  overpotential  was  found  to  decrease  with 
increasing  temperature.  At  a  higher  operating  temperature,  the 
electrodes  became  more  reactive,  and  thus,  the  activation  over¬ 
potential  decreased  and  more  current  is  generated.  This  is  related 
to  using  the  temperature  dependent  exchange  current  density  in 
the  anodic  Butler-Volmer  equation  for  evaluation  of  the  activation 
overpotential. 

As  expected,  it  was  found  that  the  ohmic  overpotential 
decreased  drastically  with  increasing  temperature  from  0.1823  V  at 
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Fig.  7.  Distribution  of  temperature  in  the  DCFC  (operating  cell  voltage  VceU  =  0.4  V). 


873  I<  to  0.1343  at  973  I<  because  the  high  ionic  conductivity  of  the 
electrolyte  was  very  sensitive  to  temperature.  Therefore,  in  order  to 
improve  the  performance  of  an  intermediate  temperature  DCFC 
operating  at  873-973  K,  the  electrolyte  should  have  a  high  ionic 
conductivity,  and  the  thickness  should  be  kept  smaller  than  800  pm 
to  reduce  the  ohmic  overpotential. 

It  is  well  known  that  the  DCFC  performance  depends  greatly  on 
the  cathode  oxygen  gas  concentration  [16],  The  result  of  the  present 
model  confirms  this  dependence  as  shown  in  Fig.  11  where  three  O2 
concentrations  are  tested.  The  DCFC  performance  is  enhanced  and 
the  limiting  current  density  increases  when  higher  O2  concentra¬ 
tion  is  added  to  the  cathode  gas  since  more  reacting  gaseous  spe¬ 
cies  are  reaching  the  TPB  which  affect  the  reaction  rate  and 
decreases  the  concentration  overpotential.  The  DCFC  seems  to  have 
a  high  limit  operating  cathode  gas  composition  since  the  amelio¬ 
ration  of  performance  between  2.060  mol  m  3  and  4.015  mol  m-3 
are  limited.  Based  on  the  present  model  results,  the  optimum 
cathode  gas  composition  needed  to  achieve  high  performance  of 
the  composite  electrolyte  DCFC  system  is  CO2/O2  =  2/1  mol%.  These 
results  agree  well  with  those  reported  earlier  by  experiments  [16], 

Fig.  11  exhibits  a  multiplicity  of  current— voltage  curve  profile  at 
different  oxygen  inlet  concentration.  For  small  oxygen  concentra¬ 
tions  (0.062  mol  m-3),  a  non-standard  behavior  was  observed  at  a 
current  density  of  about  100  mA  cm~2  relative  to  the  occurrence  of 
a  non-linear  operating  behavior.  The  /—  V  curve  exhibits  a  branch 
with  a  negative  differential  resistance  (NDR)  which  could  be  related 
to  the  nonlinear  rate  laws  of  the  electrochemical  reactions  taking 
place  at  the  electrodes  or  to  the  transport  of  charged  and  uncharged 
species  within  the  various  DCFC  components  which  follows 
nonlinear  kinetics  laws,  multi-driving  force  transport  mechanisms 
and  even  two-phase  flows  [37], 

The  falling  of  the  polarization  curve  at  lower  O2  content  is  due  to 
mass  transfer  limitation  process  which  indicates  that  there  is  a 
minimal  O2  concentration  to  be  used  in  cathode  gas  for  better  DCFC 
functioning  (Co21im). 

Indeed,  the  low  oxygen  concentration  could  affect  the  kinetic 
mechanisms  (Eq.  (22))  which  exhibits  a  non-linear  behavior  and 
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Fig.  8.  Carbon  dioxide  concentration  profile  produced  over  the  anode  sub-domain 
(operating  cell  voltage  l/ceii  =  0.4  V). 

leads  to  this  non-standard  shape  of  the  I—V  curve.  This  observed 
non-standard  behavior  may  be  also  related  to  the  electrical  con¬ 
ductivity  of  the  electrolyte  which  depends  strongly  on  the  fuel  cell’s 
temperature.  The  heat  released  by  the  electrochemical  reaction 
increases  the  cell  temperature  and  hence  the  electrolyte  conduc¬ 
tivity  which  leads  to  a  higher  rate  of  migration  of  C02~  ions  to  the 
anode  side  and  a  higher  rate  of  carbon  oxidation.  For  a  given  cell 
voltage,  this  leads  to  a  higher  cell  current  and  to  the  generation  of 
more  heat.  This  positive  feedback  between  ions  migration  and  heat 
generation  may  cause  this  steady  state  multiplicity  in  the  DCFC 
system.  It  seems  that  there  is  presence  of  a  self-accelerating  effect 
analogous  to  water  formation  in  PEMFCs  which  is  at  the  origin  of 
the  steady-state  multiplicity  including  the  hot  spots  in  SOFCs  and 
wet  spots  in  PEMFCs  [37,38], 

The  exchange  current  density  is  concentration  dependent  as 
discussed  earlier.  CO2  concentration  is  raised  to  negative  power 
of  -1.25  and  02  concentration  is  raised  to  a  positive  power  of  0.375 
as  given  by  Eq.  (22).  Increasing  the  concentration  of  O2  increases 
the  local  cathode  current  density  and  hence  decreases  the  cathodic 
activation  polarization  as  shown  in  Fig.  11.  Prins-Jansen  et  al. 
[39,40]  have  performed  impedance  analysis  and  have  concluded 
that  the  MCFC  resistance  rises  with  C02  concentration  and  de¬ 
creases  with  a  high  O2  concentration,  which  well  agrees  with  the 
obtained  modeling  results. 

The  DCFC  shares  the  same  electrochemical  cathode  mechanism. 
Various  mechanisms  have  been  proposed  for  the  cathode  reaction 
in  MCFC.  The  oxygen  reduction  mechanism  under  the  MCFC  cath¬ 
ode  conditions  has  been  described  as  several  elemental  reactions  as 
follows  [41]: 

(48) 


Fig.  9.  Oxygen  concentration  profile  over  the  cathode  sub-domain  (operating  cell 
voltage  Vo-n  =  0.4  V). 
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Fig.  10.  Temperature  effect  on  current  density,  a)  at  973  K,  b)  -m-  at  932  K,  c)  -4- 
at  873  K. 
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Fig.  11.  Cathode  gas  composition  effect  on  the  DCFC  performance,  a)  at 
021n  =  4.157  mol  nr3,  b)  at  02i„  =  2.060  mol  rrT3,  c)  at  02in  =  0.062  mol  rrT3. 


02-  +  C02(e)<-C032-  (54) 

It  is  composed  of  two  slower  reactions  (Eqs.  (51)  and  (52))  that 
limit  the  oxygen  reduction  rate.  These  reactions  are  the  reduction  of 
the  peroxide  ions  and  the  reduction  of  the  superoxide  ions. 
Experimentally,  if  the  oxygen  reduction  mechanism  is  limited  by  a 
superoxide  or  a  peroxide  reduction  reaction,  then  one  semicircular 
arc  should  be  observed  on  the  EIS  spectra;  if  it  is  limited  by  a  mixed 
peroxide— superoxide  reduction  mechanism  with  similar  superox¬ 
ide  and  peroxide  reduction  rates,  there  should  be  two  semicircular 
arcs  more  or  less  superimposed  [40],  The  prediction  capability  of 
theses  mechanisms  is  rarely  discussed  in  literature.  Therefore,  the 
two  most  common  reaction  mechanisms,  namely  the  peroxide  [4] 
and  the  superoxide  mechanisms  [32]  are  employed  and 
compared  in  Fig.  12.  Both  mechanisms  have  the  same  general  form 
but  vary  by  the  exponents  of  the  reaction  orders  in  Eq.  (22). 

Obtained  results  reveal  that  the  superoxide  mechanism  causes 
higher  DCFC  performance.  This  is  explained  by  the  tendency  of  the 
reaction  rate  (or  volumetric  current  density)  to  increase  with 
greater  magnitude  of  the  negative  exponent  of  the  C02  concen¬ 
tration  (n  =  -0.75)  and  larger  magnitude  of  the  positive  value  of  02 
concentration  (r2  =  0.625)  in  Eq.  (22). 

Moreover,  numerical  simulations  indicate  that  at  700  °C,  the 
peroxide  mechanism  is  more  predictive  to  the  experimental  anal¬ 
ysis  than  the  superoxide  mechanisms.  So,  the  cathode  mechanism 
seems  to  be  under  peroxide  mechanisms  electrochemistry  control. 

Similar  to  electrolyte  conductivity,  the  cathode  electrochemical 
reduction  rate  vary  significantly  with  different  materials  use.  The 
state-of-art  cathode  material  in  MCFC  is  NiO  with  an  exchange 
current  density  of  8.1  A  m  2.  The  tested  composite  cathode  is  made 
of  lithiated  NiO  [17], 

Fig.  13  presents  the  DCFC  polarization  curves  for  various  cathode 
exchange  current  density  values  at  700  °C.  As  seen  from  the  plot, 
varying  j0c  has  a  significant  effect  on  the  polarization  loss.  As  joc 
decreases  the  cathodic  activation  overpotential  increases  as  a  result 
of  increased  kinetic  resistance.  Consequently,  the  performance 
drops  drastically  to  69  mA  cm-2  at  lower  exchange  current  density. 
Therefore,  the  development  of  scientific  research  on  cathode  ma¬ 
terials  properties  enhancement  needs  to  be  more  intensive. 

As  shown  in  Fig.  14,  an  increase  in  the  anode  specific  area  in¬ 
creases  the  DCFC  performance  at  an  operating  temperature  of 
700  °C.  The  rise  in  the  specific  surface  area  enlarges  the  anodic 
reaction  site  and  consequently  permits  an  efficient  wetting  ability 


of  carbon  particle  with  molten  carbonate.  This  leads  to  a  faster 
carbon  oxidation  reaction.  Thus,  the  energy  barrier  needed  for  the 
occurrence  of  the  reaction  named  activation  polarization  is 
decreased  resulting  in  improved  cell  efficiency. 

This  parameter  varies  dependently  on  the  carbon  fuel  sample 
thus  it  is  very  important  to  describe  the  anodic  active  reaction  sites. 
The  activation  polarization  at  1/th  -  Weil  =  0.92  V  decreases  from 
0.43  V  at  60  cm  1  to  0.23  V  at  600  cm4.  This  is  again  consistent 
with  the  findings  from  Fig.  6  indicating  that  the  anode  activation 
polarization  had  the  most  significant  contribution  to  the  overall 
performance. 

On  the  light  of  the  obtained  results,  we  can  say  that  for 
improving  the  DCFC  performance  it  is  necessary  to  use  a  carbon 
with  a  high  specific  surface  area  to  enhance  its  wet  ability  by  the 
molten  carbonate  electrolyte. 

The  effect  of  the  electrolyte  porosity  on  the  cell  voltage  curve  is 
shown  in  Fig.  15.  The  use  of  low  electrolyte  porosity  has  a  significant 
impact  on  fuel  cell  performance  by  reducing  the  overall  current 
limit  without  showing  the  mass  transfer  limitation  falling.  Using 
different  electrolyte  matrix  porosity,  the  mass  diffusion  limitation 
process  does  not  appear  in  the  DCFC  polarization  curve  due  to  the 
fact  that  in  DCFC  concentration  losses  are  negligible  compared  to 
the  activation  and  ohmic  losses.  By  increasing  the  electrolyte 
porosity,  the  diffusion  of  carbonate  ions  (CO2-)  to  the  anode  reac¬ 
tion  site  becomes  easier.  So,  the  anodic  reaction  becomes  more 
active  and  then  reduces  the  anodic  activation  losses.  At  higher 
electrolyte  porosity,  the  concentration  overpotentials  became  clear. 
As  can  be  seen  in  Fig.  15,  the  higher  the  electrolyte  porosity  the 
greater  the  power  delivered  by  the  DCFC.  We  can  notice  that  the 
performance  of  the  DCFC  seems  to  be  effectively  improved  by 
adjusting  the  microstructure  of  the  composite  electrolyte.  This 
adjustment  could  be  made  through  trying  several  preparation 
methods  of  the  composite  electrolyte  including  Sol-Gel  process, 
Glycerine— Nitrate  combustion  process,  carbonate  co-precipitation 
or  solid  state  reaction.  The  review  of  the  influence  of  the  above 
cited  preparation  method  on  the  final  composite  electrolyte  pow¬ 
der  is  reviewed  by  several  authors  [42,43], 

The  DCFC  performances  increases  as  the  cathode  inlet  gas  flow 
increases.  Hence  the  polarization  decreases  with  the  rise  in  the  gas 
flow  velocity  as  shown  in  Fig.  16.  The  cell  performance  is  enhanced 
by  the  using  larger  amount  of  reactant  because  more  ionic 
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Current  density  (mA  cm  2) 

Fig.  13.  Comparison  of  the  DCFC  performance  for  varying  cathode  independent  tem¬ 
perature  exchange  current  density,  a)  -0-  at  7  A  m~2,  b)  at  50  A  irT2,  c)  -m-  at 
500  A  nr2. 


generation  in  the  TPB  can  be  reached  via  accelerated  cathodic 
reduction  and  anodic  carbon  oxidation  electrochemical  reaction 
rates.  The  high  velocity  reduces  the  concentration  gradient  be¬ 
tween  bulk  and  superficial  layer  which  decreases  the  concentration 
polarization  effect. 

5.  Conclusion 

In  this  study,  a  2-D  numerical  model  of  a  direct  carbon  fuel  cell  is 
developed  to  simulate  the  cell  performance  by  taking  into  account 
of  the  electrochemical  reaction  mechanisms,  mass  and  heat 
transfer  and  electrode  materials  characteristics  and  processes.  This 
proposed  model  allows  the  prediction  of  the  effect  of  the  operating 
conditions  (temperature,  cathode  kinetic,  matrix  porosity  and 
anode  specific  area)  on  the  DCFC  performance.  It  permits  also  the 
investigation  of  temperature  and  species  concentrations  profiles 
within  the  DCFC  based  on  composite  electrolyte.  Experimental 
validation  of  the  DCFC  performance  fueled  by  a  mixture  of  graphite 
and  molten  carbonate  is  also  performed  showing  a  good  agreement 
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Fig.  15.  Variation  of  the  DCFC  performance  with  the  electrolyte  matrix  porosity,  a)  -#- 
at  30%,  b)  -m-  at  20%,  c)  at  10%. 


with  an  absolute  average  deviation  of  about  9%.  Results  show  that 
the  linear  tendency  of  the  polarization  curve  is  not  solely  related  to 
the  ohmic  loss  control  but  also  it  is  a  result  of  a  combined  effect  of 
the  anode  and  cathode  activation  loss.  The  temperature  distribu¬ 
tion  through  the  DCFC  computational  domain  varied  by  15  °C 
among  all  the  cell  sub-domains  due  to  the  near  zero  entropy 
change.  Comparison  between  the  cathode  kinetic  mechanisms  in¬ 
dicates  that  the  peroxide  mechanism  predicts  the  DCFC  polariza¬ 
tion  curve  accurately.  Results  show  that  the  high  cathode  inlet  gas 
flow  velocity  is  a  key  parameter  accelerating  the  electrochemical 
reaction  in  the  cathode  TPB.  Simulations  are  conducted  to  evaluate 
the  effect  of  the  anode  specific  surface  area.  It  is  found  that 
increasing  this  parameter  can  enhance  the  cell  performance.  The 
produced  current  density  is  also  significantly  improved  at  a  more 
porous  electrolyte  matrix.  Further  work  is  needed  to  improve  the 
capabilities  of  the  present  model.  In  fact,  we  project  to  include  ef¬ 
fects  of  the  non  local  thermal  equilibrium  within  the  porous  media 
(between  the  solid  and  gaseous  phase).  This  model  will  also  be 
extended  to  simulate  three-dimensional  effects  in  order  to  capture 
geometrical  phenomena  effects. 


